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A DTSGUSSICN rp DYNA^^IG STABILITY MSAS'REMSNT TECHNIQUES 

V.’ayne Johnson* 

Ames Research Center and 
U.S. Army Air Mobility R&D Laboratory 


SUMMARY 

Techniques for the neasurement of the dynamic stability of linear 
systems are discussed . Particular attention is given to an analysis of the 
errors in the procedures, and to methods for calculating the system damping 
from the data. The techniques discussed include: transient decay, moving 

block analysi'-, spectral analysis, random decrement signatures, transfer 
function analysis, and parameter identification methods. The special 
oroblems of rctorcraft dyna.mic stability testing are discussed. 


INTMnXJCTTC'N 

The flutter testing of an airplane or helicopter requires a method 
for measuring the dynamic stability of the aercelastic system. The objective 
of the test is to establish the sxability level throughout the operating 
range of the aircraft, "he capability to accurately determine the stability 
trends is necessary to safely conduct such a test. Thus a method is required 
to reliably retermine the frequency and damping of all important modes of 
the system, using measurements of the system response to existing or prescribed 
inputs. Lynanlc stability measurement is a particularly difficult task for 
rotorcraft, where many '’agrees of freedom are involved and the process and 
measurement noise ’eve Is can be very high. 

A number of techniques have been implemented or suggested for measuring 
dynamic stability. These techniques may be usefully classified as transient, 
no-input, input-output, and parameter identification methods. Transient methods 
use the decaying transient response of the system to determine the dynamic 
stability. The moving block analysis is a particular nrocediore for analyzing 
transient data. The no-input methods use the response due to existing random 
disturbances of the system (such as aerodynamic turbulence), with averaging of 
the response to reduce the errors in the estimates. Examples of no-input methods 

♦Research Scientist, Large Scale Aerodynamics Branch. 
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are spectral analysis, correlation, and randon decrement signatures. The 
input-output methods use the response due to a measureable external input to 
determine the system characte-^istics alone, as well as tc reduce the error. 
Examples of input-output methods are cross- correlation and traasfer function 
analyses. Parameter identification metho’s use statistical estimaxion 
procedures to obtain a best estimate of the system characteristics from the 
measured response, .ill of these method’s have been used, in rotorcraft or 
airplane flutter testing, with varying degrees of success depending on the 
aircraft involved and the test environment. 

The stability measurement procedure must perform with some level of 
process and m.easurement noise present, .i primary concern is how such noise 
reduces the accuracy of the measurements, hence introduces error into the 
determination of the system damping. This report present.^ an error analysis 
of the various stability measurement methods. Some of the results are new 
(such as the random decrement signature error analysis and the discussion of 
parameter identification techniques for rotorcraft), while others are available 
in the standard literature. A unified treatment of the error analysis allows 
a consistent comparison of the accuracy of the available stability measurement 
techniques. Thus this rejxirt provides a guide for choosing the technique 
appropriate to a particular experiment. 

The basic stability measurement procedixre in all cases gives only an 
estimate of the response, in either the time or the frequency domain. Further 
processing of the data is required to determine the system parameters from the 
response, in particular to obtain the damping ratio which gives the quantitative 
level of stability. Thus for each procedure some alternatives for data processing 
are discussed. The following stability measurement techniciues are studied in 
detail: transient decay, spectral analysis, random decrement signatures, 

transfer function analysis, and parameter identification methods. First however, 
a mathematical definition is .“riven for the aeroelastic system considered. 

LINEAR SYSTEM DYNAMICS 

To examine the dynamic stability measurement techniques, a linear 
system is con.oidered, with excitation by various control inputs and external 
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disturbances. The response of the system is measured, perhaps with significant 
measurement noise such as vibration due to rotor imbalance or vibration from 
the engine. The excitation of the system can be urovided by existing unknowr. 
disturbances such as aerodynamic turbulence, or by a measureable external 
input applied snecifically to determine the dynamic stability. The system 
motion is therefore described by linear, time-invariant differential equations, 
of the form: 

X = Fx + Gu 

y == Hx + V 

where x is the state vector and y the observation. The vector v is random 
measurement noise; and u is a random input exciting the system, either an 
existing disturbance or an external input. The matrices F, G, and H are 
constant since the system is time-invariant. The solution for the response 
to excitation by u, with initial conditions at t^, is 


f (i.-to') 


V ^ Cr u 

Wo 


1^(0= H X it) 4- V 


The stability of the system is determined by the eigenvalues of F . The 
eigenvalues usually occur in complex conjugate pairs, of the form 

X I - g ^ ' 


where Ov.. is the natural frequency and ^ the damping ratio of the mode. 

The mode is stable if ^ > 0, and is an exponentially decaying oscillation 
for 0< ^ < 1. It is usually the lew darii-d nodes of the system which are 
of interest, i.e. 0 <5< 0.1 . Note that if -A_ is the diagonal matrix 
of the eigenvalues of F , and II the matrix whose columns are the corresponding 
eigenvectors, then ~ For further information on the dynamics of 

linear systems, the rearer is directed to references 1-3* 


Regardless of the stability measurement pioced’jre used, a basic 
requirement is that the system be controllable and observable. This means 
that the input to the system must sufficiently excite the modes of interest, 
and the modes must be observable in the response of the variables measured. 
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TRANSIENT DECAY 

The transient notion of a dynamic system is composed of exponentially 
decaying oscillations of each mode. A particular measurement is often 
dominated by one, low-damped mode. From the oscillation period and the 
decay rate of the trace, tlie frequency and damping of the mode may be 
estimated. The typical procedure involves establishing a large amplitude 
sinusoidal motion by means of an external excitation at the natural 
frequency of the particular mode of interest. The excitation is stopped, 
and then the subsequent transient motion is analyzed (often by hand) to 
determine the damping. 


Consider the linear, time-invariant system 

X = Fx + Gu 
y = Hx + V 

where x is the state and y the measurement; u is process noise such as 

aerodynamic turbulence, and v is measurement noise. It is assumed that 

u and V are random disturbances with zero mean. The motion following 

time t , where the initial conditions x(t ) are established by the external 
o ^ o' 

excitation, is 

^ Gru &‘c 




Hx -♦-V 


The first term in x is the transient we wish to observe. The second term 

is the noise in the trace due to the disturbances occuring after t . 

o 

The influence of this noise is the primary concern with this method . 


Statistics 


The expected value of the observation is just the transient 


response ! 




The external excitation must establish a large, non-zero initial condition 
for the modes of interest. The variance of the observation is 
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where it is assumed that u is white noise with Eu(’t^ )u (f^) = QS(ti-Tj), 
so also u(t>t ) is uncorrelated with x(t ). The measurement noise 

O Ip ^ 

correlation is Ev(’C,)v (Xj) = R(x,-l£). It has been assumed that the 
initial conditions x(t^) are deterministic, not random, variables. 


Considering a single degree of freedom system, so the integrations 
are easily performed, the expected value and variance of the observation 
are : 



JU. X, 


V(X--to') 

9 . 



-*t 

St. St -V <- 




t-') 




2-V 


-V r 


So the normalized variance is 


1 

€. 





«" V — 2X iX-Xo^ 

A.' 



Note that the <^teady sta .e variance of the rtate x wiih only the excitation u 

= g^q/(-2 X). so 




/V — N -2V(-t“-ko') V 


Xo' 


t 2 2 

At t = t^ the only contribution is from the measurement noise, (. = r/h x^. 
As t increases, the noise due to the ciisturbance u becomes more important. 


The transient is decaying, while the motion due to u quickly builds up 
to the level V_ 
exponentially : 


to the level V^. Thus the ratio of the noise to the transient grows 


, -2X(-t-XO ^ 

■ \ u - o 

X* * 
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To determine the damping ratio it is necessary to observe the transient 
while it decays to a fraction f of the initial value (f = .3 to .5 typically). 
'I'he maximum error which must be dealt with is thus 



So C is of the or'’er of 2 to 3 times the ratio of the rms response 

due to u and the initial response amplitude. This result holds for the 

case of more than one degree of freedom also. Note that the steady state var-iance. 

' ^f the resoome of the system to the f'isturbance u alone is given by 

- ( rx + 

from wh LCh it follows that V X. 

y 

To reduce the error in this procedure, one can only increase the 

amplitude of the initial conditions x(t^), or reduce the noise soiirces 

u and V, The measurement noise v is frequently small, or at least under 

some control by the experimenter (e.g. the rotor l/rev vibration can be 

reduced by improving the track). The process noise can also sometimes 

be controlled, e.g. by conducting a flight test under calm atmospheric 

conditions only. In wind tunnel experiments however, often there is little 

which can be done to reduce the turbulence significantly. Increasing the 

initial excitation may not be practical either; the largest value possible 

2 

is often used to stctrt with. In summary then, the noise parameter 

often can not be controlled by the experimenter. This parameter is typically 

large (say .5”). so that the error in the transient trace is a major problem. 

Data Processing 

The frequency and damping are often determined by hand, from an 
oscillograph t-^ace rf the response. It is assiimed that the motion i.mr 
the form of a decaying oscillation of a single mode: 

X — V o ^ * t-%.4 C w« 'Ai- 'jj* "3 

Then the frequency is given by the period of the oscillation, to = 1/T (Hz), 

To determine the damping ratio, consider two peak amplitude measurements x^^ 
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and Xg, which are n oscillations apart. Then 

*x=». ^ ^ 

*'-S‘ 

and 

4m / >1 ( 

^ “3^^ ^ 

For example, ^ equals .11 divided by the number of cycles to one-half 
amplitude. 


x^exp(- ^i^t) 


£)Jl^p /fir 

2yr to 


where cJ is the oscillation frequency in Hz. The damping ratio then nay 
be calculated from the slope of the log-decrement; the previous method if 
just a special case, where the slope is calculated from the two points 
Xj and x^. One data processing procedure is to plot In x vs. t (removing 
the mean value first if necessary). The slope of the envelope of this plot 
gives the damping ratio ^ . In general, determining the envelope of the 
transient is a difficult task to mechanize however. 


In general, the envelope of the oscillation x = 
gives In X = In x^ - I^us^t, or 

/<Qx 


5 = 


CO 


I + ( 


A different approach is to find a least-sq>jared-error fit to the 
data for the curve 

y ■= > ^ ^ ^ ^ ) 

Identifying the parameters r- and is a nonlinear problem however, so an 
iterative solution is required. There have been some successful applications 
of this approach. However, convergence problems are very likely with a 
high noise level. 
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Movliyr Block Analysis 

The moving block analysis is a particular (iata reduction technioue 
for transient notion. The Fourier transform of a block of data from t 
to ’c+T gives a function y(*c). The magnitude of the sprectrum line at 
the natural frequency w. plotted vs. *c gives then the damping ratio. 


The moving block analysis is defined by the linear operator Lj 


L = 


_L 

z-^ 



•o t. 


5)t 




C*-) 1 

-t»r 


(Sit 


So y(*t) = Ly(t). The expected value of y is 


E^= = 


Ptt-t o') 

m CL )Co 


-A-tt— -to') 

= UM Lo 


- HM 




T- y (t-» t -“to - .. 

Ot J V 


= M e. M 


r r 1 


The dependence on »c is e ^ , so y is still a measure of the transient 
motion. The last factor in Ey is indepent'ent of *t , and contains the 
only dependence on the spectrum frequency . This factor is a maximum 
when the frequency XL is equal to the natural frequency of V . When 
S~L is near the natural frequency of a mode, the q’jantity 


- I 

•L-n C V - Jri > 


I - 

■2-.T- 


is large; away from the resonance it is small. Thus the largest line 
of the spectrum "y will t.e at the natural frequency of the mode dominating 
the response y. Considering only the line at n - therefore, and 
neglecting all modes but the dominant one (with eigenvector 1 j) , we Jiave: 

-f - 
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^51 


5Lxi»^ 


=. ^ 


l-e 




z*tr^l>v, 




magnitude of this line Is: 

I ^ 

j'xis 


2:= CjO<« 


s e. 


- 5vJ^^T-tO 


:p. 


an exponentially decaying function of T . The moving block analysis 
thxjs results in the envelope of the decaying oscillation directly; this 
is the principal advantage of the technique. 

In summary, the moving block analysis involves applying the Fourier 
transform to the block of data from t to T. + T. Usually the FFT algorithm 
is used, with a sampled time series of data. The r ignitude of the line at 
the natural frequency is obtained; if the data is not too noisy, this 
line will have the largest magnitude in the entire spectrum, and the frequency 
of the mode may be determined by searching for this peedt rather than being 
input. A plot of ly\ vs. x gives then the envelope of the decaying 
transient; the slope of Injy^ vs. t gives the damping ratio ^ as 
discussed above. Since the procedure gives the envelope directly. It 
is simple to mechanize the calculation of the damping ratio, such as 
by finding the least-squared -error fit of a straight line to Injyj vs. t* 

This is a significant advantage compared tc the processing required to 
find the damping ratio from the transient motion y. 

In choosing the block size T, it can not be too small or the operator 
reduces to simply L = i/?'K ; and if ^w^T is small, then the S2. = 

line of the spectrum does not have the largest magnitude. So there must 
be a significant decay of the oscillation in the time T. A block length of 
3 or 4 periods is probably satisfactory. Tc increase the frequency resolution, 
the FFT should not be applied to the block of length T, but to one several 
times as long, filled with zeros after the data. 
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Now consider the einror involved In the moving block analysis . 

The linear operator L can have no fundamental influence on the error. 

The normalized variance of y(*t ) nust have basically the sane value as 
that of y(t), l.e. error level determined by the parameter V Working 
in the frequency domain helps, since only the noise around the frequency 

is important. This may be expec+ad to reduce the effect of the measurement 
noise significantly. The operator has little effect on ohe process noise 
however. The response of the system to u will be mainly a superposition 
of oscillations at the natural frequency , which then is transmitted 
to y along with the transient oscillation. 


The variance of y is 

V 




The evaluation of Vy is complicated by the operator L. Considering only 
the domi. nt mode in the response, the expected value and variance are 
approximately 

_ 




= He. 


The normalized error is 




L ( 


) Je, H 










which is indeed essentially the same result as for the normalized error 
of the original transient trace. The minor improvement in the error 
statistics due to the moving block operator is not the reason for its use. 
The principal advantage of the procedure is that it. extracts the frequency 
and exponential envelope from the decaying oscillation, which gi'eatly 
facilitate: the mechani'/.ation of the data processin^^ task. 
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SIBCTRAL ANALYSIS 


To improve the accuracy of the meaFurement of the svstem respor.ae, the 
data mubt he averaged in some fashion. Spectral analysis is a fundamental 
technique for measuring the random response of a system. The excitation 
of the sys ten is provided by existing, random disturbances, such as 
aerodynamic turbulence, inly the response y of the system is measured. 

Let Sy be the autos pectrum of a particular measured time series of d.ata. 

Then an estimator of the true spectrum of the response is given by the 
<^nsemble average over K calculations of the spectrum: 


S 






The spectrum S 

y 


may be obtained by either digital or analog processing. 


In general, th’ process noise exciting the system may consist 
of multiple inputs (e.g. several gust compor.:ints ) , so the true output 
spectrum is related to 'he input spectri by the expression 


where is the transfer function of the response y to the i-th input, 
and is the ^ unknown) input autcspect-^um . lor the linear system 
under considerat’on, H, is a rational function, and the denominator 
is the same for 'ill inputs (it is the characteristic equation). Thus 


s. = ^ 




He are primarily interested in the poles oi the characteristic equation, 
which give the damping cf the ry- tem. The output autos])ectrum is 
a direct nearuire of I). Ibere will be a resonant peak of S near the 

i 

ne.tural frequency oJ a low-damped mode, h‘- long as there is no correspending 
zero of Njj near the pole. Such a numerator rero may arise if the mode 
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is not observable in y, or if the input spectrum is varying greatly in 
the vicinity of the resonance. Therefore the mode must contribute 
significantly to the measured response y, as usual; and it is necessary 
that the input spectrum be reasonably flat near the resonance. The latter 
restriction is not too severe since the resonances of low-damped modes 
are very sharp. 

upectiTal analysis and correlation analysis are entirely equivalent, 
the former in the frequency domain and the latter in the time domain. 

Thus the correlation techniques will not be discussed here; their statistics 
are similar to those of spectral analysis. With a digital processor 
and hardware FFT units, the spectral analysis is more convenient; it is 
also better suited for data reduction, such as finding the damping ratio. 


Statistics 

A 

An analysis of the statistics of the estimator of the autospectrum 
is available in the literature (see references 4 and 5)* 1’he expected 
value and variance of S are 

sm 

V sC*-') 

< 

where is the frequency resolution in the spectrum, and K is the number 

of average i. These results are essentially independent of the statistical 
properties of the input and output, hence are applicable to the response 
to general random disturbances. The normalized bias error is 

, s'" 

The most critical case is at a resonant oeak, where the greatest accuracy 
is req red, and S'' is largest. Assumiing that in the vicinity of the peak 
the spectrum may be approximated by that of a single degree of fi'eedom 
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system, one obtains S/S*' = i 5**^' • ^ given bias error, the 

following frequency resolution is required : 


Then j gives = .02, which is essentially negligible 

bias error (twice that resolution is probably satisfactory). Note that 
the half-power bandwidth of the peak is = 2.50^ , so this resolution 

corresponds to covering the bandwidth of the peak with 5 spectral lines. 


The bias error is easily made negligible by a proper choice of 
the fr . ency resolution. The remaining error is the normalized 
variance 1 _ . a _ l 


= [cirs^ V “ >nr 


The error is inversely proportional to the square-root of the number 
of averages, the standard result for sample means. The estimate of the 
spectrum may thus be made as accurate as desired, by increasing the 
number of averages. 


An important parameter of the spectral analysis technique is the 
total record length required: 

_ ■■ ^ — 

So the total time over which the data is collected must be increased as 
the damping, frequency, or error is decreased. This is a fundamental 
result for the amount of data required to define the system from the 
output due to rand urn dis curbances , applicable to all no- input techniques. 

Data Processing 

In using the cuipu.,' -■‘ctrum to determine the system properties, 
it is generally assumed that the in] ut spectrum is flat in the vicinity of 
the resonances. Since the low-damped modes are of most interest, the 
assumption is only required over the narrow frequency range 2^ . 
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Assuming that the spectrum in the vicinity of a resonance may be approximated 
by that of a single, second order modes 


S 



CAa* S’VowaaT 


it follows that the damping ratio may be determined from the spectrum 
peak and an integral through the resonant peaks 
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RAN1X)M DECREMENT SIGNATURES 
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The method of random decrement signatures (randomdec) is a procedure 
developed by Cile (reference 6) for analyzing the response ol a linear 
system to random disturbances. It is designed to give an estimate of 
the impulse response of the system. The randomdec estimator "y is 
defined as the ensemble average of the transient response to existing 


random disturbances: 


— 


_ -L £ 


L±'> 


V- » 


The key to the prpceaure is to select each ensemble yj^ with the same 
initial conditions by an appropriate triggering method. Then all records 
have an identical transient due to the initial conditions, while the subsequent 
noise averages out. 


Consider a linear system excited by the a«kn*wn random disturbance 
u; the response measurement is y, and the measurement noise is v. 

The system is described by the differential equations; 

X = F'x + Gu 
y = Hx + V 

The solution using the initial conditions at time t^ is 

xct'^= V 4- 

^ to 

= H xCt') -V V 


The observation y(t) is a random process, depending on the random variables 

u, V, and x(t ). The expected value of the randcmdec signature is 
0 



■yT ^ 


P(t-0 

t Ev^ 




/-t-to) 


X o 
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where it has been assumed that u and v have zero mean, and that the 
initial conditions are stationary random variables with expected values 
= Ex(t^). Thus y is an unbiased estimator of the system transient 
response. 


In order for a particular mode to be analyzed, it must be observable 
in y, and excitable by u. The initial conditions x(t^) are due to the 
random disturbance u over the tine up to t^, so the frequency content of 
u must be sufficient to excite the modes of interest. This is a fundamental 
requirement of methods using existing random disturbances of the system, 
arising here through the initial conditions at time t^ . Furthermore, in 
order for a particular mode to be present in y, it is necessary that the 
initial condition for that mode have a non-zero average value over all 
the ensembles, i.e, Bx(t^)^0. The response x(t) is a random process 
with mean 

Ex = F"x(Su) = 0 

So it is not possible to arbitrarily select the ensembles. The result 
would be and hence E7 -■ 0. It is necessary to trigger on the response 

in such a manner that Ex(t^) f 0 for the modes of interest. Note that 
triggering on y=0 is not possible either, since it gives zero initial conditions 
(it would be possible to trigger on y=0, with positive slope only, and 
hence obtain a non-zero average value of the initial velocity; the variance 
of the initial velocity would be large howe/er, leading to a large error 
in y as discussed below). The usucil procedure is to trigger the start of 
the ensemble on a fixed, non-zero level of the response y, with either 
positive or negative sloi>e. 


As usual, the ergodic hypothesis is used, and the random decrement 
signature is obtained by averaging over records which are sequential in 
time : i K 






It is assumed that the records are uncorrelated, which is generally 
satisfied if there is no overlap (ref. 6). With digital sampling and 
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processing of the data, the requirement that each record be triggered on 
the basis of the response means that the records must be non-overlapping. 
Using records with some overlap would reduce the total sample time some 
(a factor of 2 is about the maximum possible reduction), but Increases 
the correlation of the records. 


Statistics 


The variance of the random decrement signature y is 


= E 5^ - (Ef^r 








The ensembles are uncorrelated, except for the initial conditions. The 

initial condition x(t ) is a random variable with mean Ex(t ) = x and 

2 ^ o' o 

variance Vx^= Ex(t^)x (t^) - , determined by triggering on the response 

y in some fashion. Note that the response of the system to the random 

disturbance u has zero mean and variance = Ex(t)x (t), which are not the 

same as x^ and Vx^ due to the action of the procedure for selecting the 

initial time t . 

o 


Now assume that u and v are uncorrelated with each other. Let 

T ' T 

3v(*c, )v (Tj) = ; and assume white noise for u, so Eu(tj)u (tr^) = 

Then also u(t > is uncorrelated with x(t^). Then the 

variance of 7 is 







, ,T- 




e 




e. 


<9t vT 



The term in brackets is the same error as was found for the transient 
decay trace. It is the response of the system to the random disturbance 
u subsequent to the time t^ where the initial conditions are determined. 
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This term is zero at t = (except for the measurement noise contribution), 
and increases with time until it reaches the steady state level given by V^. 
The relative size of this term increases exponentially as the transient signal 
decays. The total normalized error due to this source is approximately 


C = 


C£5V 


VC 




So the ensemble averaging reduces this error, inversely proportional to 
the square-root of the number of averages. This error source can be 
made as small as required by Increasing the number of averages. 


The remaining term in Vy depends upon the variance of the initial 

2 

conditions, Both this term and (Ey) decrease with time at the same 

rate, so the normalized errcn is always about the same value: 



Thus the normalized error of y is equal to the normalized error in the 

Initial conditions x(t ). The averaging process doss not reduce this error 

0 

term, so the accuracy of the signatui'e depends upon keeoing the variance 
of x(t^) as low as possible. 

The variance Vxq depends on the procedure for triggering the start 
of a sample record , Zero error = 0 requires that the same initial 
conditions x^ be produced exactly by every triggering event. An important 
observation concerns the influence of the matrix H in Vy. Since (at t = t ) 

Vy = , it follows that the variance in the initial conditions of the 

unobserved states do not contribute to Vy'. The triggering process, based 
on the measurement y, can not be expected to produce a specified initial 
condition for an unobservable state, for such states the initial condition 
will have a zero mean, and a variance equal to the value in V^. This 
large variance does not contribute to Vy however, so it is of no concern in 
the randomdec procedure. An important example is the velocity of the 
degrees of freedom of an aeroelastic system. The state variables consist 
of the displacement and velocity of the degrees of freedom, but generally only 
the displacements are measured. The velocities are unobservable states 
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then; the velocity initial conciitions will have large variance when the 
triggering process uses the measurer^ response y, but only the errors in 
the displacement initial coniiltions contribute to Vy, 

Consider a single degree-of-freedom , second order system, with 
states X and x. The response x is measured, and the saunple record is started 
by triggering on the level of x. Then Vy is just the variance of the initial 
condition on x alone. An efficient triggering mrocess should give a small 
'■Xq for the displacement. Therefore this example should have low error 
in the random decrement signature y due to the initial condition errors. 

Obviously the triggering procedure Is a central element in the 

calculation of the random decrement signature. A procedure is required 

to establish the initial conditions with the following properties. The 

expected value Ex(t^) = x^ mast be large enough for the important modes 

that they have signific?.nt excitation in the measured transient. The 

variance of the initial conditions \'y /x ' must be small for the observable 

^0 o 

states; this is essential for low error in the signature. Finally, 

2 

either V must be small, or a la‘»"ge number of averages must be made 
so that the error due to the random disturbances after t is small. ( f 
course can not be very small, since the sane random input u produces 

both ' and x(t), of which x( t ) is a member. In sunnarv then, x should’ 
be large both to reduce the errors and to Increase the magnitude of y. 

If too large a value of is chosen for the triggering process however, 
the correlation of the records increases. Cole (ref. 6) suggests, on the 
basis of an analysis of the statistics of the method, tl^t the triggering 
use a level y = 1.0 to 1.6^ (where is the ms value of the rertonse). 

ftinimizing the error due to the initial conditions variance requires 

that the triggering process select t on the basis of the measurement y 

c 

such that the respoase x has the prescribed value. There are two basic 
reasons why this may be difficult tc accomplish, both of which must always 
be present to some extent: the existence of measurement noise, and systems 

with more than one mode obser'/able in the measurement. 'Measurement noise 
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prt:;ent in y will give a lower bounc.' on the variance \' . ' nliJce the ^/K 
term, this effect of the measurement noise does not average cut. Depending 
on the syrtem, this may be a serious noise soiirce or s negligible one. Note 
Uvi t errors in the triggering process h;ive an effect similar to the measurement 
noise . 


A usur'il case involves a single measurement of the response of a 
multi-mode system. If more than one no^ft is obs.ervable in y, then obvic’isly 
it is not possible to uniquely detenaine the values of the obrei’vable states 
and thus trigger the record solely base'- on the single value y(t). In 
general the initial confections can not be detf-innined if the dimension of 
y is less thiat the dimension of the observable states in x^. In practice, 
only tliose states which contribute -.ignif icantly to y are of concern, 
and f req iently there- i s only one such state contributing to the measurement. 

An example of where more than one mode is obse7~vable is the neasurenent of 
the tip motion of a cantilever bean (e.g. a wing or rotor blaf’e). In such 
a case the frequency range may be restricted, or only a single more excited, 
so that there is no practical difficulty; or in fact it nay be necessary to 
determine the initial conditions of all the nodes based on tha one measurement. 
M, is (X)sslble to obtain an estimate of the state xf t) from limited noisy 
meas'trements y, if use is ma'le of the past history of y. This is a classical 
filterin;; problem; the solution in this case i;- the Kalman filter, which 
deperifis on the characteristics of the system being investigated. It is 
possible to Identify the r;ystem as the measurement?', proceed, and thus update 
the inter model. This becomes .a rather involved procedure simply to 
trigger the p.-birt of the record. 

In summary, it in often ]x>srlble to trigger on y with low variance 
of the initial conditions of all obsrr’/ahle states, hence with low error 
in the ranriomdec slgaature jr. V/ith a high order system, or with measurement 
noise the triggering tank becomes more difficult, and the error in ” 
increases. In the extreme case, the triggering procedure becomes the 
[)r'oblem of identifying the state x( t) from noisy measurements of an unknown 
linear process -- essentially a definition of the entire dynamic stability 
measurement ta:.k. 
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It r.houlf’ be noted that with -rmltlple stater observable in y, all 
are required to take the prescribed initial conditions at the tri(^gerln{7 
instant. Vhe tire interval between such occurrences v.iil increase as the 
number of states increases, hence the tot-al testing time (not the sampling 
time) will increase. Alternatively, fer a given testing time the error 
must be allovied to increase, either by relaxing the trigger criterion an-' 

•'0 increasing or by baking fewer averages. 

Vne above result for ’'y v;as based on white noise for the r^in'or. 

disturbance u. 'Mie result will be similar for a general snectrur. , 

)u (*Cj) “ T, - *C j ) , which indeed is more realistic than white 

noise. The basic influence of averaging in reducing the noise is unchanged , 

although the number of averages required for a specific voice of error 

will be influenced somewhat by the spectr'im of u. If the input spectrim 

is low at the frequency o"^ the node to he investigated , it will be necessary 

to increase the number of averages to conpensote for the lew signal to 

noise ratio there. An imnor+ant consideration for toe random decrersr.t 

procedure is the effect of handwidth-linite ’ ^nput disturbances on the 

initial conditions. The random d ist'urbance u is the source of x( t) , f'^om 

which the initia] conditions must be chosen. Hence the frequency content 

of the input has a direct influence on 3 x( t ) and . 

o' 

I'inally, consider the sarnie time require', Cole gives an e''tinate 
of T based on the stxatistics of y". ’ or the problem of near ■»'ing the 

(hamping ratio, the result is 

where £ is the normalised error of the (iamping estimate. This is the 
same result as for spectral analysis (altho;:gh the i n ter pret~ tier/' of 
£ are somewhat different in the two cases). The amount of ''ata required . r. 
extract information about the system witn a given accuracy is a fundanenbal 
characteristic, independent of the procedure used to extract the infr ■'■natior.. 
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As with spectral .••nalysis, for u given T it is necessary to trar’e 
off the number of averages and the frequency resolution, '^he frequency 
resolution of the randondec signature Is rtetermineo by how many perlo'^r 
of the oscillation ire measured. If there are n cycler with frequency 
In the rit/naL’.re, -then the lengt.h of each -rcor'^ is t = ri/t.>vv • 

Vhe P' mber of aveoge'" ' v then 



- ith ^ - .01 to .02, 6 = .1 to .2, and n = 4, the n'jnbei" of averages 

required is I - 1000 to SOCO. 'Jole suggests the ui-e of n = h, and a sample 
rate of r ~ (so the signature conrist^s of samples) . Tynical 

applications of the ran-’or '’ecreuent proce'^ure requi'^e K - 4000 to "000 
averages for satisfactory accuracy preferences 6 an' ?). 

Data i recess i rig 

'.he random d;ecrenent signature y is. an estimate of the system transient 
response. Thus the proco' ures riiscusse'’ above for eterminlng the den' ing 
ratio from the decaying transient trace are applicable, vith ■'mpreve''' 
accuracy due to the averaging. eterninlng the 'anrinr an' l>equency 
by means, of a least-rqua red-error fit cf a 'decaying oscillation to the data 
is more successful than with the unaverciged trac , hue to the reduction in 
the error of the estimate of the transient. jh.ang (’-ef, 7) devel-iped such 
a data Processing technique for the ran'ondec netho'. 

In conclusion, the spectral analysis, correlation, an-'* rani-’or! 
decrement techniques are all menrerr of the general class of methods 
for analyzing the proper‘‘ies of linear dynamic systems utilizing existing 
(unknown) randora ('isturb.7tnces, , with averaging tri reduce the e^Tor. 

As such, all of these techniques have basically the same statistical 
characteristics. The choice of method depen'^'r on the expc-rlneat, the 
data processing equipment, an' as. often on the preferences of the 
experimenter. Correlation is the classical technique, .v'ith current 
'digital processors, and F.-'T hardware, it is more useful and convenient to 
work In the frequency domain with spectral analysis. The rindom decrement 
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ctatir.ticr. are r.omewhat rpecial in their '^epen’ence on the variance of 

the initial con<'n.tlonG , anr’ hence on the trif^frerinf; I'roceso. The cennitivity 

of Vv to meacurenent noine, an'^ the (difficulties cf the inulti-mot^e 

care arc of ^reatect concern. The spectr'.n of the output f^eDenio -'irectly 

on ■frequency content of the in'out, while the rcin'onHec signature 

’eucndr. on the input opectrun ln'’irectly , particularly throu/;h x an'* ' 

o 

All procedures are eroentlally in'^epen(’ent of the level of the input, as 
expected for a linear system. The nomalization is in'olicit In the 
rancdon’ec trif^gerinf, an'’ 3 may he normalize’ by '’jviding by ^ ”’o . 

All cf the procedures require further procer.ring to extract t'-e damring 
ratio and other para-neterr. In this regar' the soeclrf-l analysis ir mere 
convenient thfin correlation or inndcndec tec'^'niaues , which give the response 
in the time domain. 
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•nAN:JfST FUN ITK U AN'ALYGIS 

Stability analyris techniques which rely on existinc, unknown 
TP-nrion dirturhancee to excite the systen have the baric limitation that 
the measured response repen’s oii the properties of both the system an^' 
the unknown input, ’iy measuring, the input as well as the responre, it 
is poFsible to determine the characteristics of the system alone, as well 
as reduce the error in the analysis, .,'sually it is necessary to apply 
a measurable external innut speci fieri lly to excite the system. The 
cross-spec tra between the input and output, in addition to the autospectra 
are found. The ratio gives the transfer function H (O) between the input 
w and the resTX>nre y, from which the system jnrameters may be determined. 


Consider now a linear system wit i i epu'’ and output near 'ur erne nts : 

X = F X -r Gu 3w 
y - X ^ v^ 

7 . - W d- V^ 

The actual system response and input are x and w; y and z are the measured 
response and input, Inclu'ing ran'’on measurement noir e. .he vector u 
is the process nolr.c, an unJ-'nown random dirtu’^bance of the system. The 
transfer function between x and w is 

U(u'i = Q = oiw-py'r 


Let S be the cross-snectr'im between m, and y based on the k-th record 

(kl Ik) 

of data, and and 3^ the cor rcrs ponding antes jiecta. '"hen an emtimaioir 
• r’ the transfer function H(u) is given by the ratio of the averaged 
v^ross-spectrijn and input autospectrun : 
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The arvilyi',iE of the statlftlcr of the er-tlmator H of the 
Ky^+en frequency respcrue ir available in the literature (see references 
/J’ anr* 5)» fot' the spectral analysis technique the variance of the 

A A ' 

soe''tr=‘ estinators- S and S wene r’iscussed, iho statistic? of their 
ratio \i are different hov.’cver. reflecting the 'act that not .lust the input 
or output is being nearured, but the correlativon between the two. Bias 

A 

errorr- in H a*'e usioally such lens than the random errors (as long as the 
frequency resolution is chosen properly), so they are not discussed, 

A 

The analysis of the variance, rt H gives the following result. 


The probability is oc t'lat the magnitu-^e of the difference 

A 

between the estimate H(tj , and the true value H(t<i) is greater than 
the function r(ui). Thus. v;ith confidence '^1 - cy ) H lies within a circle 

A 

of radius r about H on the plane of InH vs. Be’.i. i.'o 


1 < r 


where ^ ^ 

^ -P- 

A 

Here L is the nun'oer of averages: is the estimate of the coherence 

function between the input and output 


A 

y 






and i* is the /alue of the :• distr 
freedom at the probability level 
Note that ideally |h\‘' - .1 /ll : 

* y 7 . 

^2 5 a fu- O 


Ibuticn witg !" and Jl(K-l) d 
' \ - o<> ) , i .c . i rob{ ■ > 
am* for larre i' fin’ 

~ — 


egreer of 
) - cx ). 


Hence the normalize'^ error of tl.e e: +imate of t'^e transfer f.-nctlon la 




For a confidence level of 05 :' ( o< = . 05 ), find F ■= 3*1» a-t the 97.5 

anc’ 99 levels, F ^ 3*^ and 4,8 respectively. Assuming a Gaussian distribution, 

the variance (one standard deviation) corresponds tc a confioenra level 

A 

of 68 , which gives F = I.I 6 . Hence except for cbs effect of , the 

X A. / " 

normalized error is € = 1/K, which is the same result as for the variance 
of the spectrum alone. 


A 

The coherence function is a measure of how well z and y are 

correlated, '.'ith no noise in the system, would be exactly equal to 1. 

z. 

3o (1 - is. a measure of the effect of noise oa the calculation of 

the response due to a partictilar input. For the system described above, 
assume u, v. , aa’ v^ are uncorrelated random noise; let H be the frequency 
raspoase of x due to w, and the frequency response of x due to the 
d i s turba nee u , The n 
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The first tern is the ratio of the system resTjor.se due tc :he disturbance u, 
to the total response; it is always ler.s than 1 since the lat+- is the 
sum of the response to u and w. This, tern nay he •^educed b’ incn'j<sinr 
the level of the measure-'' input w relativ-" to the unknown inp...u u. The 
last two terms are the ratio f the msasv.rer.ent noire to signal for th- 
ou tput and input. 


It is concluded that for given mimber of averages (hence flxn ' 
sample time), the -=rror in the transfer fu^'ction is substr^ntially less 
than the error in the individual spectra, "'he error nay be direc* y 
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reduced by increasing the level of the measured excitation w. The improved 

error statistics compared to spectral analysis or equivalent techniques 

u 

is a particularly important factor for full-scale testing, where the 
frequencies are low and tha« the sample times required tend to be long 
(recall that T = 2/C'^-C*')* The other principal advantage of the 
transfer function analysis is that it produces more information about 
the system than does the spectral analysis of the response alone. The 
result is independent of the input spectixan (always assuming that the 

input sufficiently excites the relevant modes), 

n 

IData Processing 

Processing of the frequency response data to obtain the damping 
is similar in some respects to the processing for spectral analysis, although 
of co’-urse the assxanption that the input spectrm is flat in the vicinity 
of the resonance is not required, A typical procedure involves approximating 
H near the resonance of a low-damped mode by the frequency response of a 
single, second-order degree of freedom: 



^ - w*- i 


Tt follows that the daunping ratio may be calculated from integrals of the 
measured transfer function about the resonance; 
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The greatest effect of other modes on the resonance is usually a phase 
shift. The assumed single degree- of- freedom response has a phase Z.H = -90^ 
at the natural frequency, but the phase of the measured H can be shifted 
significantly by nearby nodes. Thus the measured frequency response phase 
•Ti Id be shifted as follows before applying the above expression to 
calculate T : , ^ 

^ u. = u ^ (-1 ^ 

where Hp is the measured response at the resonant peak. 








An alternative approach is to use the methods of Kennedy and 
Pancu (reference 8). Cn the ImH vs, ReH plane, the transfer function 
describes a circle in the vicinity of each resonance. Near the natural 
frequency, the ideal response of the single degree-of-freedom ^system is 
approximately 
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The radius of this circle thus gives the damping of the mode. The 
procedure nay be implemented as follows. First the data should be 
shifted to phase -90° at the peak frequency, as indicated above. Then 
the circle 

+ z, - 1= 

on the H plane is fitted to the measured transfer function, using the 
least-squared -error method to calculate a and b. Cnly a few points around . 
the peak should be used, say 3 to 5 either side* and a weighting 
function should emphasize the errors for the points nearest the peak. 

Then the damping coefficient is 




To find the damping ratio however, it is necessary to further consider 
the variation of the phase with O at the resonance. A least-squared -error 
fit of the line 


to the data gives the slope c, and hence the damping ratio ^ “ 2/c. 

A simpler procedure results if the inverse of the transfer function is 
considered. The ideal, single degree-of-freedom response H ^ is exactly 
a parabola on the complex plane: 

-1 -1 2 

Thus the curve (ReH ) + a(lmH ) = b, may be fitted to the data, using 

the least-squared -error approach. Then the damping ratio is given by 
5 = (4ab)"2. 
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Another approach is to fit the data to a high order rational 
function, e.g. by least-squared -error methods. The roots of the denominator 
polynomial then give the stability of all modes of the system obrervable 
in the response. The assumed representation of H depends nonlinearly 
on the parameters however, so an iterative solution is required. Convergence 
problems are possible with such a procedure, especially if there is 
significant noise in the measured transfer function. 
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Parameter identlficatlor techniques for dynamic stability measurement 
require further development. It is possible to define the problem for 
rotorcraft testing however. An algorithm is required to estimate the 
frequency and damping of all important modes of the system from measurements 
of the various input and output signals. For rotorcraft there are usually 
many degrees-of-freedom in the system, and a high noise level in the 
measurements. It is necessary to identify the structure an'^ st‘».tistics 
of the noise as well as the system for good accuracy. The assumptions 
about the system must be minimal. It is seldom possible to make assumptions, 
about the model without compromising the effectiveness in identifying 
unforeseen stability characteristics. In particular the order of the system 
m'ist be identified, i.e. the number of significant nodes involved in the 
response. The algorithm must give an accurate estimate of the damping of 
the system, with assured convergen''e and computational efficiency. It must 
operate on-line, providing the estimates in a reasonable time after the 
collection of the data for all input combinations, with the calculations 
performed by a limited capacity computer. 

It is assumed that the system is described by linear, time- invariant 
differential equations, excited by various control inputs and aerodynamic 
turbulence. The turbulence is defined as the response of a linear process 
excited by white noise. The system response and input are measiored . The 
measurements can have three kinds of noise: bias error (balance or zero 

Ehifx,), random error (not necessarily white noise), and discrete frequency 
error (correlated with rotor azimuth, e.g. l/rev or w/rev vibration). The 
input signals may often be considered noise-free. Thus the following 
description of the system may be used. Modal coordinates are used 
since the system is then defined in minimal form, and directly in terms 
of the eigenvalues X (the eigenvalues are probably complex however, sc 
actually a real form of the following description is required). The 
model is; 
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system: q = 

output measurement: y = 

input measurement: z = 

turbulence model: u = 

Here y and z are measurements of the response and the input, with noise z 

V (bias, random, and discrete components). The aerodynajnic turbulence 
is u, excited by the white noise v^. 

The data consists of the measurements of the response and the 
input, probably in the form of sampled time series. Sequential measurements 
are made for the various input combinations (including no input). Then the 
parameter identification algorithm is applied. It is necessary to identify 
the system (-A^ , G, E, H, order), the turbulence model (D, order, v^), and 
the measurement noise model (v^, v^, v^). The eigenvalues _A- are of 
primary concern of course. Clearly the magnitude of the parameter identification 
task is very great. 

There is a serious difficulty with measurements in the rotating 
frame. With independent blades (i.e. shaft-fixed rotor dynamics) the 
measurements and the equations defining the system can be in the rotating 
frame. With coupled motion of the blades and the rotor support, the states 
and equations must be defined in the nonrotating however, including those 
for the rotor (otherwise the system is described by periodic coefficient 
equations, needlessly in this case). For example, the rotor flapping notion 
should be described by the coning and tip-jiath-plane tilt motions of the 
rotor rather than by the flapping of the individual blades. If only 
measirrements in the nonrotating frame are made, there is no difficulty, 
except that due to limited information. If measurements are available 
for the motion of all the blades of the rotor, then they can be transformed 
by analog or digital processing to the nonro a.ting degrees-of-freedom. 

But very seldom are rotors instrumented for the stability measurement 
requirement, so almost always just one blade will be instrumented. It 
would be very valuable to have an algorithm capable of using the single 
blade measurements in determining the dynamic stability of the system. 


■ -f ' ; 


A. q + Gu + Ew 
Hq + Voj + Vo, + voj,^ 
''in ’ID 
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The methods which have been discussed so far (spectral analysis, and 
transfer function techniques) can not use single blade measurements, 
except for the case of independent blade dynamics. These techniques 
require measurements representing the response of a time -invariant system, 
so they must have measurements from all the blades or they are not able 
to use the additional information. Parameter identification techniques 
however have the potential for using the limited information in the single 
blaf^e measurements. It would not be as accurate as having measurements 
of the motion of all the blades, but it would be better than ignoring the 
rotor motion entirely. 


Consider the revisions of the measurement representation required 
to account for iiaving only rotating measurements for a single blade. 

Assume for now that there is only one measurement y of the blade motion. 

Let yj^ be the corresponding N measuren ents of the states in the nonrotating 
frame (" is the nunbc" of ’-r'+or bla'^es). So y^ is re?ated to the modal 
coordinates by y.^p = Kq as usual. Jiow y is related to the nonrotating states 
by a linear (time-varying) transformation y = h^^y.^; i.e. the response y 

All i 

is due to the response of all the components. For example, the flap 
motion of the m-th blade is given by 

' * V» S I V, 

-t- peC-tl 

where =ri,t d- 2'rr/N, and ^n>;z is present o»Ji^ N is even. Thus 

h is the 1-by-N vector: 

I 

C-OS 

(-O'* J 


It follows that y is related to the modal resTXjnre q by the expression 

T 

y = h Hq + V 

o 

This result is easily extended to the general case. Let the first n_ 

r 
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elements of y be the meacurements in the fixer’ frame, and the last n^^ 
elements of y be the measurements on the rotating blade. 'Ilien 

y = THq + 

where 

~r rz 

The matrix h and therefore T are known. The only effect of using the 
rotating blade measurements is an increase in the dimension of the matrix 
H which must be identified. Note that H is just the matrix which must 
be found if all N of the nonrotating rotor states are measured, rather 
than the rotating degrees-of-freedom of the single blade. ITie matrix 
T (i.e. h ) accounts for the reduction in measurements (and infomation) 
when only the single blade instrumentation is used. 

For rotorcraft dynamic stability measurements, it is also necessary 
to consider the periodic coefficient case. There are two usual sources of 
periodic coefficients; the aerodynamics of the rotor in forward flight, 
i.e. 0 ( .5 in practice); or a two-bladed rotor on an anisotroptic 

support. In the periodic coefficient case the system is described by the 
same model as above. The eigenvalues are given by -A- as us’jal. Now 
however the matrices K, G, and S are periodic functions of time (the period 
is 2vr/Nr2. if in the nonrotating frame, and ?.tT/rL in the rotating frame, 
where iT. is the rotor speed). The parameter identification problem is 
increased therefore, since there matrices must be found over the entire 
period now. Undoubtedly the required sample time will be increased. 

It would probably be desirable to synchronize the sampling of the data 
with the rotor speed (with the trim rotor speed, since there are cases 
where the dynamic response involves significant rotor speed Tjerturbations) . 
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aNCLUDTNG RRHAnKS 

A nunber of methor's for measuring the dynamic stability of linear 
systems have been discussed. Particular attention has been given to an 
analysis of the errors involved. It was found that a significant 
reduction in the error is achieved when averaging is introduced, and again 
when measureable external inputs are used. The results are a guide to 
the design of ex]>eriments and the choice of the measurement technique. 

No universal conclusion about which technique should be used is poss.ible. 

The choice depends on the experiment, what noise is iiresent, what error is 
acceptable, and on what data processing equipment is available, -'-■z' 
rotorcraft testing the dynamic stability measurement is a critical factor 
however, so generally the most acc'orate method possible is desired. 

The state-of-the-art of dynamic stability testing in the aircraft 
inr’ustry has been described in a number of surveys (see for example 
references 9-13) • The problem involved is usually the flutter testing cf 
the aircraft wing or tail in flight. Some form of f-^equency-response 
analysis procedure is common. There are many variations, but a typical 
procedure involves swert-sine excitation, '’igital analysis of the response 
(analog analysis is also still common), with the damning det mined from 
the circle on the H plane (see the above discussion, and ref. 3). 

The techniques discussed may be implemented in many ways. Among 
the significant options are: the response may be in the time or freauency 

domain; the input may be swept-s.ine or ranc’om excitation; analog or 'igital 
pj'ocessing may be used; and the data reduction may be accomplisher manually, 
with a special pur'oore processor, a mini-conputer , or with a large capacity 
computer. The discussions here have dealt with those ciiaracteristicr. 
which ore fundamental to the problem of dynamic stability measurement, an^ 
hence are independent of the method of implementing the procedure'-. 

finally, this report has been limited to an examination of the basic 
characteristics of the methods for measuring dynamic stability. ■ urther 
development of techniques for rotorcraft testing requires computer simulations, 
and additional applications to stability measurements in rotor experiments. 
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